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1. Introduction 

The ultimate goal of the ultra-relativistic heavy-ion collisions being currently carried out 
at CERN is the observation of the phase transition from hadronic matter to the so-called 
quark-gluon plasma [1, 2]. One expects that during such a deconfinement phase transition, in 
which the structures of individual hadrons are destroyed and replaced by a highly-energetic 
system of quarks and gluons, chiral symmetry is additionally restored. Lattice simulations of 
QCD [3] indicate that the chiral symmetry restoration phase transition and deconfinement 
phase transitions may indeed occur at the same temperature. Looking for the signatures of 
these two phase transitions is one of the most challenging problems of high-energy nuclear 
physics, both from the experimental and the theoretical point of view [4]. 

Since the matter produced in high energy nuclear collisions lives only for a short while, 
it is natural to expect that its space-time evolution proceeds far away from equilibrium. 
Consequently, there exists a growing interest in applying and developing transport theories, 
which are suitable for the description of non-equilibrium processes. A kinetic theory based 
directly on QCD has been formulated in papers by Heinz [5], and by Elze, Gyulassy and 
Vasak [6]. However this approach, with a few exceptions [7], is still far removed from practical 
applications. 

In the context of the two phase transitions that have been mentioned, it is desirable 
and important to be able to include some of their aspects into the kinetic approach, which 
can otherwise be based on some phenomenological assumptions. The deconfinement phase 
transition seems at present to be too difficult to deal with, since we still do not yet understand 
this phenomenon fully. On the other hand, the chiral phase transition has been studied in 
terms of effective chiral models of QCD, for which phenomenological Lagrangians can be 
given. Such an approach has been quite successful in clarifying this transition. 

One such effective approach is based on the Nambu - Jona-Lasinio (NJL) model; within 
its framework we can study the temperature (density) dependence of various physical quan- 
tities like, e.g., the quark condensate, the pion decay constant, or the quark and meson 
masses [8, 9]. It is also possible to study the change of the thermodynamic properties of the 
quark-meson plasma as the system passes the chiral phase transition [10]. With few excep- 
tions however [11, 12, 13], the considerations based on the NJL model are usually restricted 
to equilibrium situations. We refer the reader to reviews on the NJL model given in [14]. 

In the present paper, we wish to investigate how the concept of chiral symmetry can 
be explicitly included into the quantum and classical transport theory based on the NJL 
Lagrangian. In our investigations, we shall restrict ourselves to the mean field or Hartree 
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approximation, since the interesting phenomena (like, e.g., spontanous symmetry breaking) 
already take place at this level. Our approach is based on the spinor decomposition of the 
Wigner function and follows the treatment of [15] and [16]. However, in contrast to [15] we 
do not neglect the spin degrees of freedom and take into account all the components of the 
spinor decomposition. On the other hand, we differ from [16], where the spin dynamics is 
fully discussed, in that we study a system that is governed by a different type of interaction. 

A transport theory for the NJL model has been already formulated by Zhang and 
Wilets [11], where the Keldysh closed-time-path formalism combined with the effective- 
action method has been used to derive the kinetic equations. A certain class of solutions 
of the classical transport equations found in [11] has been studied in both [12] and [13]. 
Nevertheless, Ref. [11] does not discuss how the symmetry of the underlying theory (i.e., 
the chiral symmetry of the NJL Lagrangian) should reflect itself in the transport equations. 

In this paper, in trying to clarify this point, we derive the general form of the transport 
equations (in the mean field approximation) which are chirally invariant. We discuss the 
properties of these equations and then study their classical limit by performing an expansion 
in h for all the functions appearing in the spinor decomposition. Throughout our approach, 
we investigate and check the consistency of our equations under the assumptions of chiral 
symmetry. In particular, this implies that the axial current should always be conserved. For 
the spin evolution equation, as will be seen, this turns out to be a non-trivial requirement. 

Our paper is organized as follows: In Section 2 we define the NJL model in the mean 
field approximation. The Wigner function and its spinor decomposition are introduced in 
Section 3. In Section 4 we present the chirally invariant quantum kinetic equations. Their 
classical limit is discussed thoroughly in Sections 5-8. In Section 9 we study the case when 
chiral symmetry is explicitly broken by the presence of nonzero current quark masses. We 
summarize and conclude in Section 10. 

2. The NJL model in the mean field approximation 

In our study, we consider the following Lagrangian of the NJL type, 

£ = j, (i@ - m ) V> + G + (^TbVO 2 ] , (1) 

where ij) is the Dirac field, m is the current quark mass, and G is the coupling constant. 
For simplicity, flavor and color degrees of freedom have been neglected. In the case m = 0, 
the Lagrangian (1) is invariant under £7^(1) transformations. This symmetry is, however, 
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spontaneously broken by the ground state. In a more realistic approach to quark matter, one 
should consider the Lagrange density which is invariant under SUa{%) transformations, with 
Ua{X) explicitly broken due to the presence of anomalies and with the spontanous symmetry 
breaking of SUa{%) [14]. Nevertheless, we expect that the basic problems concerning the 
formulation of the chirally invariant transport theory can be already studied in the framework 
of a simpler theory, like this based on (1). In what follows we shall always assume thatm = 
except for Section 9 where the consequences of setting m ^ are taken into consideration. 

By introducing the fields a and 7r defined as 



*=-2Grl>rl>, ir = -2Grl>i-y 6 il>, ( 2 ) 
we can recast Lagrangian (1) into the form 



C = ipi flip — aipip — ■Kipi / y 5 ip — — . (3) 

4Gr 

The equivalence of (1) and (3) becomes clear if one checks that the variation of (1) with 
respect to if; gives the same equation of motion for ij) as the variation of (3) with respect to 
ij),<T and 7r, namely 

[i p- a(x) - i<y 5 ic(x)]ip(x) = 0. (4) 

In the following, we shall restrict ourselves to the mean field approximation, i.e., the operators 
<x and 7r will be replaced by their mean values 



<x — > a = (a) = Tr(/5<r), 7r — > 7r = (7r) = Tr(/J7r). (5) 

where p is the density operator and Tr denotes the trace over all physical states of the 
system. The particular form of p is not important for our considerations. We only note that 
in typical situations, p is either an operator projecting on the parity invariant ground state or 
it describes parity invariant statistical ensembles. In these two cases (7r) = 0. Nevertheless, 
in more general situations (e.g., when the system is out of equilibrium) we cannot exclude 
the case (7r) ^ 0. 

The Lagrangian (3) is invariant under the £7y(l) transformations, if; — > ij)' = exp(— 
which leads to the conservation of the vector current, 

d»V»(x) = 0, V(x) = {i,(x)^i,(x)). (6) 
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In addition, as was stated before, it is invariant under £7^(1) chiral transformations 

V> -> V>' = exp(-»7 B ^)V>, (7) 
<x — > cr' = <x cos % — 7r sin %, (8) 

7T — > 7r' = 7T COS % + CT sin %. (9) 

Consequently, the axial current is conserved as well, 

d^x) = 0, A»(x) = {j>(xW' l6 1>(x)). (10) 

In general, for a non-equilibrium transport theory, one requires two independent Green 
functions to describe the system. In the mean field approximation, however, one is sufficient. 
We thus choose the Green function 

GZp{x,y) = {MvW.{x)) (11) 

to be the fundamental quantity in our description. One can easily notice that it satisfies the 
same equation as the field ij}{x) does, i.e., in the mean field approximation we have 

[ifl- a ( x ) - ; 757 r(z)] G<(x,y) = 0. (12) 

The mean fields cr(a;) and tt(x) are determined by the Green function G < (x,y) through the 
relations 

a(x) = -2Gtr G < (x,x), tt(s) = - 2Gtr »7 B G !< (s, s), (13) 
where the trace tr is taken over the spinor indices. 

3. Wigner function and its spinor decomposition 

For the description of non-uniform systems, it is convenient to introduce Wigner trans- 
forms. This allows one to make an easier physical interpretation of the components of the 
Green function and facilitates performing the classical limit. In the case of the Green func- 
tion G < (x,y), its Wigner transform is obtained by introducing the center-of-mass coordinate 
X = [x + 2/)/2, the relative coordinate u = x — y, and by taking the Fourier transform with 
respect to u. Following the notation of [17], we define 

W a p(X,p) = J d\eir- u G<p (X + ^J- |) . (14) 
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In Eq. (14) we have explicitly incorporated the Planck constant h, since later on we wish to 
investigate the classical approximation. 

The Wigner transforms of the derivative of a two-point function, df(x,y)/dx ,J ', and of 
the product of a one-point function with the two-point one, f(x)g(x,y), are given by the 
expressions 

dJ ^ L - (-^ + ^)/(X, P ), (15) 

f{x)g{x,y) - f(X)g(X,p)-jdJ(X)dZg(X,p), (16) 

with the notation = d/dX^,d£ = d/dp^, and where in the second relation (16), only the 
derivatives of first order in d p and dx are kept. This approximation will always be used in 
the following. 

Using Eq. (12) and the properties of the Wigner transform (15) and (16), we arrive at 
the following equation for the Wigner function 

^ + jd^j ltl - a(x) + ycU(x)d; - ; 75 7r(x) - l lh d.-K{x)d; w(x, P ) = o. 

(17) 

Eq. (17) together with Eqs. (13) and (14) form the system of the coupled equations that we 
are going to study in detail in this paper. We observe that under the chiral transformation 
(7), the Wigner function changes according to the prescription 

W _> w> = exp(_i 7B |) W exp(-t 7B |). (18) 

Hence, one can check that Eq. (17) is invariant under transformations (8), (9) and (18) 
and thus that the first order in the derivative expansion of the Green function is symmetry 
preserving. 

The Wigner function W a p(X, p) is a 4 by 4 matrix in spinor indices and can be expanded 
in terms of 16 independent generators of the Clifford algebra. The conventional basis consists 
of 1, 275, 7^, 7^75, and which we denote in what follows by I 1 ;. In this basis, the Wigner 
function has the form 

W = T + i-y 5 V + + 7^75^ + ^<V- ( 19 ) 
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The properties 7°^/ = W and 7 rt 7 = V i (the first one results from (14), whereas the 
second one is a simple consequence of the definition of T;) indicate that the coefficients in the 
spinor decomposition !F(X, p), V(X, p), V M (X, p), A^(X, p) and S^ v (X,p) are real functions. 
We note that the spinor decomposition technique has been already used in the formulation of 
the transport theory for QED [16, 18, 19] and for QHD (so-called quantum hadrodynamics) 
[15]. 

The vector and axial currents, see Eqs. (6) and (10), are related to the functions V M (X, p) 
and A^(X, p) by the simple relations 

d A p 



V*(X)=4f—^V''(X,p) (20) 



and 



A»{X) = -4/ ^LA^p). (21) 



Therefore V M (X, p) and A^(X, p) can be interpreted as the momentum densities of these two 
currents. For the physical interpretation of the other components we refer the reader to Refs. 
[16, 18]. Under the chiral transformation (18) the coefficients of the spinor decomposition 
(19) change as follows 



T - 


■* r = 


J~ cos x + V sin %, 


V - 


■* v = 


—J 7 sin x + V cos Xi 


- 


■* v = 


v„ 


A.^i 


-> A' - 


A 


Sfj,V - 


■* S' = 

flu 


= S^v cos x + «$U sin X> 


Sfj,V - 




= -S^ sin x + cos x 



(22) 



where S**" is the dual tensor to S**" , namely 

1 
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4. Kinetic equations 

Substituting expression (19) into Eq. (17) gives the following system of the coupled 
equations for the coefficients of the decomposition (19) 





-aF + TvV = 


-y (d„ad;r - d„ird;v) , 


(24) 


- iK»Ar 


- aV - tv J 7 = 


-j (d v ad;v + d^d;?) , 


(25) 


K^T + %K V S VI1 - 


aV^ + iivA^ = 




(26) 


iK^V - K V S^ - 


aA^ + iirV 1 = 


-j(d v ad;A>*-id v Tvd;v>*), 


(27) 


i(K»V v - K^^-e^ K a A p - 


-irS^ + aS^ = 




(28) 



In abbreviation we have introduced here the notation X" = p" + yd". Eqs. (24) - (28) 
should be supplemented by the formulas that determine the mean fields. They follow from 
(13) and have the form 

and 

Using the transformation properties (8), (9) and (22) we can check that this system 
of equations is chirally invariant. Strictly speaking, each of Eqs. (24), (25) and (28) is 
separately a chirally invariant equation. On the other hand, Eqs. (26) and (27), after a 
chiral transformation, form linear combinations of each other. These combinations imply, 
however, that each of their components must vanish separately. In consequence, we can 
treat the system of the two equations (26) and (27) as a chirally invariant expression. In the 
analogous way, one can check that the system of Eqs. (29) and (30) is chirally invariant. 

For further analysis, it is convenient to discuss separately the real and imaginary parts 
of Eqs. (24) - (28). The real parts give 



(29) 
(30) 



+ = 0, 



(31) 



^d^-aV-irF = 0, (32) 

p^-^S^-aV^ = -^d v Kd;A^ (33) 

-^V-pJ^-aA 11 = -\d v -Kd v v y\ (34) 

-\{^ V -d'V^-e^PnAp + aS* 1 " -ttS^ = 0. (35) 

One can notice that integrating of Eq. (32) over the momentum gives 



hj-0^d^(X,p) = O, (36) 

where the definitions (29) and (30) have been used. This equation is nothing other than the 
statement that the axial current should be conserved (10). We thus see that after making 
the gradient expansion, this conservation law is still included in the transport equations. 
The imaginary parts of Eqs. (24) - (28) yield 



^v, = -\(d v *d;r -d^o;v), (37) 
= ^ (d v ad;v + d v Kd;T) , (38) 

^F + p'S^ + tvA^ = -\d v ad v p V^ (39) 
^ - 1^5"" + ttV = -\d v ad;A^ (40) 
( P »V V -p'V^-^e^dnAp = \ [d^S^ - d^d^) . (41) 

Although we have neglected the higher order gradients, Eqs. (31) - (41) are still quantum 
kinetic equations. In order to obtain the classical equations one makes an expansion in powers 
of T in h 



T = F {0) + &F ( i) + h 2 F {2) + ... , (42) 

and similarly of V, V M , «4 M , S^", tt and a. Expansions of the form (42) are inserted into Eqs. 
(31) - (41) and the expressions at the appropriate powers of % are compared. The detailed 
description of this procedure will be the subject of the next Sections. 
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5. Constraint equations in the leading order of h 



Substituting expansions of the type (42) into Eqs. (31) - (35), we find to leading (zeroth) 
order of % 



P^ - ^(0)^(0) + 7r (0 )7>(o) = 0, (43) 

o-(o)^(o) + ^(0)^(0) = °> ( 44 ) 

P^o) ~ <r W Vfa = 0, (45) 

SjfijP* - "W-Ato) = °> (46) 

-e^p^ + awSft-irSft = 0. (47) 



Correspondingly, the zeroth order of Eqs. (37) - (41) has the form 





= 


rs<$ + 7r(°U(°) 


= 




= 


p'Vlo) - P'Vfo 


= 



(48) 
(49) 
(50) 
(51) 

These are only four equations, since Eq. (37) is already of first order in %. From Eqs. (44) 
and (45) one finds 

P(o) = -Ho)^ 1 (52) 

and 

Then Eqs. (50) and (51) are automatically fulfilled. Moreover, Eqs. (43) - (45) lead to the 
mass-shell constraint for the function J r ( )(X, p), namely 



[p 2 - M 2 (X)]F [0) (X,p) = 0, M 2 (X) = a 2 [0) (X) + tt ( 2 0) (X). (54) 
From Eq. (47), we find the following expression for the spin tensor 

C**» _ _^(0) r „ _ v .p 1 Z(0) uvaf) a(o) , . 

(°) ~ M 2 \? (0) M 2 13 ' ^ ' 
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and the dual spin tensor 



Substituting now Eq. (56) into Eq. (46), and using (48) one obtains a mass-shell constraint 
for A( j(X,p) also, i.e., 



[p 2 -M 2 (X)]^ o) (X,p) = 0. (57) 

Finally, using (55), (48) and (57) we find that (49) is satisfied. In summary, Eqs. (43) - (51) 
lead to the expressions for V(o) and in terms of ^(0)5 and for and in terms of 
-Atj^y Furthermore, they lead to the two mass-shell constraints. 

The mean fields appearing in Eqs. (43) - (47), (49) and (50) are required to be calculated 
in a self-consistent way from Eqs. (29) and (30). Doing so, we find 

a {0) (X) = -8CJ ^JP (0) (X,p), (58) 

7r ( o)(X) = 8G J j^V {0) (X,p). (59) 

Due to the relation (44), one finds that these two equations are not independent and can 
be reduced to a single equation which determines the invariant mass M(X). To see this 
relation more clearly, let us define the distribution function F(X,p) via the expression 

f(x,?) = = (60) 

and the angle $(X) via the relations 



7r ( o) = M(X) sin $(X), o- ( o) = M(X) cos $(X). (61) 

F(X,p) is a chirally invariant quantity, and knowledge of it determines the value of M(X). 
In addition, it allows us to calculate the vector current density V^. However, Eqs. (58) and 
(59) do not determine separately the fields 7T( ) and 0"(o) 5 i.e., the angle $(X). This fact is 
in agreement with the requirements of the chiral symmetry of the problem: under the chiral 
transformations $(X) — > $'(X) = $(X) + % and, consequently, the absolute value of $(X) 
has no physical significance. 
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With this, we conclude the discussion of the expressions obtained to zeroth order in 

h. The derived equations all display the chirally invariant form. Let us turn now to the 
discussion of the equations which follow from Eqs. (37) - (41), considered to first order in %. 

6. Kinetic equation for the quark distribution functions 

Employing Eqs. (52) and (53) in Eq. (37), one immediately finds 

P»d»F{X,p) + M(X)d„M(X)d^F(X,p) = 0, (62) 

where we have used the definition (60). This is again a chirally invariant equation. Due to 
the mass-shell condition (54), we can express F(X,p) as the sum of the quark and antiquark 
distribution functions / + (X, p) and /~(X, p) defined by 

2E P (X) J K 2E P (X) [J y ' W J j ^ > 

Here E P (X) = y'p 2 + M 2 (X), and p is a three-momentum. Substituting expression (63) 
into Eq. (62), and integrating over p° gives 

prdJ^X, p) + M(X)d„M(X)d^f ± (X, p) = 0. (64) 

Eq. (64) displays a typical form for fermion mean field theories. In our case, however, a 
characteristic feature of this equation is the appearance of the chirally invariant mass M(X). 
We emphasize that the four-momentum p^ occuring in (64) is to be taken on the mass shell, 

i. e., in this case p° = E P (X). At the same time the distribution functions are functions of 
the space-time coordinate X and the three-momentum p with the derivative with respect to 
the energy, p°, defined to be zero. We have formally included it however, in order to display 
the relativistic invariance of the kinetic equation (64) explicitly. 

Using Eqs. (58), (60) and (63) we find that the NJL gap equation takes the form 

* {0) (X) = 4G* [0) (X) j ^jS^Jjq I 1 - / + ( X >P) - /"(*>P)] . (65) 

which determines the mass M(X) appearing in E P (X) in terms of the distribution functions 
/ ± (X, p). One can easily notice that using Eq. (59) instead of (58) we would obtain the 
same equation as (65), with the appearance of 7T( ) instead of cr(o)- 



F(X,p) = 2tt 
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Using now Eqs. (20) and (63), we find the following expression for the baryon current 



v *w = 2 lw?im t /+(x ' p) - r{x ' v) + J l - (66) 

where p^ is again on the mass shell, i.e., p° = E P (X). The constant appearing at the end of 
the square bracket can be neglected, since it does not contribute to V(X), and for V°(X) it 
introduces a constant (although infinite) charge. 

We note that Eqs. (64) and (65) form a closed system of equations: the distribution 
functions f ± (X,p) determine the invariant mass M(X) through Eq. (65), whereas the space 
dependence of M(X) determines in (64) the time evolution of f ± (X,p). In the case 7T( ) = 
this system of equations was first derived by Zhang and Wilets [11], whose calculations were 
based on the path integral formalism. Our approach represents an alternative derivation 
and generalizes their results to the case 7T( ) ^ 0. Some solutions of Eqs. (64) and (65) have 
been found in [12] and [13]. 

7. Spin evolution 

In the classical limit the spin dynamics is described by the behavior of the function 
A^(X, p) [16, 18]. In order to derive the kinetic equation for A^(X, p), we use Eq. (40) 
and substitute into it the expression for V M obtained from Eq. (33). Such a procedure gives, 
to the first order in h, the following equation 

P^d v A( 0) - * {0 )dJ$ - noA S {o) + Md v Md v p A[ 0) = 0. (67) 

Substituting the expression for the spin tensor (55) and the dual spin tensor (56) into Eq. 
(67) leads to the equation determining the space-time evolution of A^(X, p), namely 

p»d v A[ 0) + Md v Md;A[ 0) + ^ [pMfa - p v A[^ - e^d v ^p a A^ = 0. (68) 

We note that Eq. (68) is again a chirally invariant equation, since both the mass M(X) 
and the gradient of the angle $(X) are chirally invariant quantities. Due to the condition 
(48), only three out of four equations in (68) are independent. In fact, one can check that 
multiplication of (68) by p^ gives zero. 

The mass M(X) appearing in (68) has to be calculated from the system of equations (64) 
and (65). Hence, it can be treated in (68) as an externally prescribed function. On the other 
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hand, the gradient d^(X) is not known, and until we specify how to calculate it, Eq. (68) 
cannot be used to determine the time evolution of A^(X, p). Moreover, one has to check 
whether the solutions of (68) satisfy the requirement of axial current conservation (10). The 
last two points will be discussed in more detail in the next Section, where we examine the 
consistency of our all equations up to the first order in %. 

In analogy to the QED calculations [6] we introduce the spin up and spin down phase- 
space densities 

F ±S (X, P ) = F(X, P ) ± S,(X,p) ^ {x) P) , (69) 
where S^(X,p) is defined by 

S»( x ,p) = - -r. (70) 

[-A?\x,p)Al 0) (X,p)Y 

Due to the condition (48), the four- vector A^(X, p) is space-like, A^(X, p)A^(X, p) < 0, 
therefore S^(X,p)S^(X,p) = —1. The quantity S^(X,p) describes the mean spin orientation 
of the classical particles placed at the space-time point X and having the four-momentum 
p. On the other hand, the magnitude of A^(X, p), through Eq. (69), determines how many 
spin up and spin down particles are present in the system. Using Eq. (68), we find 

p v d v S 11 + MdvMd^" + ^~V il S v - e^ al3 d v §p a Sp = 0, (71) 

and consequently 



p»d„F ±s (X,p) + Md„Md;F ±s (X,p) = 0. (72) 

Hence, in a similar way to transport equations for QED [16], we find that the spin components 
decouple. Since both functions F(X,p) and A^(X, p) are on the mass-shell, we can write 
the analogous decomposition as in (63), namely 



F ±s (X,p) = 2. | 2Ep{x) f ±s (X, p) + 2Ep{x) [f Ts (X, -p) - Ij | . (73) 
Substituting Eq. (73) into Eq. (72), and integrating over p° gives 
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P»dJl(X,p) + Md^Md^fl (X,p) = 0. 



(74) 



By a direct comparison of (63) and (73), we find f ± (X, p) = | [/^(X, p) + f^ s (X, p)j . 
Therefore the functions / ± (X, p) simply represent the spin-averaged densities. This has 
been given in agreement with the normalization of (65) and (66). 

8. Consistency of the classical transport equations 

The kinetic equations for the quark densities (62) and for the spin (68) follow from Eqs. 

(37) and (40) considered to the first order in %. We have yet to examine (i) whether equations 

(38) , (39) and (41) are consistent with these to the same order and (ii) whether the additional 
information for the angle $(X) appearing in (68) can be obtained. We note that this kind 
of checking has been done in the case of QED in [16]. 

First, in order to illustrate that there is no contradiction between (37) and (39) we 
substitute the quantities A^ and S VI1 obtained from Eqs. (34) and (35), respectively into Eq. 

(39) . In this way, we find the following equation 




Since we are interested in retaining terms to the first order in % we can use the relations 
(44), (45), (54), and the kinetic equation (62) to find that (75), and consequently (39), is 
fulfilled. 

In an analogous fashion, using the kinetic equation for spin (68), the constraints (48) and 
(57), and the definition (61), we can show that Eq. (41) is fulfilled again to the first order in 
h. To do this, we multiply both sides of (41) by e a p^. Using the formula for V" from (33), 
we find that the multiplication of the left-hand-side of (41) gives the expression 

- ^ {e aM » [a 7 5JJ - d lHo) d;A\ 0) ] - m [d a A^ - dpAp) } . (76) 
On the other hand, the multiplication of the right-hand-side of (41) yields 
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ft [d^wd;^ + d lHo) d;s^] . (77) 

Using now the expressions for S^J and Sj®), after a rather lengthy calculation in which Eqs. 
(48), (57), (61) and (68) must be used, we find that the two expressions above are equal to 
each other. Consequently, Eq. (41) is fulfilled up to the first order in % as well. 

In this way, we have shown that Eqs. (37), (39) - (41) are satisfied in the first order. 
To prove this fact, we have used Eqs. (31) - (35), treating them as if they were satisfied up 
to the first order. Is this procedure admissable? The answer to this question is yes. Only 
T and A* are independent functions (with A* restricted by the axial current conservation 
(36)). Eqs. (32), (33) and (35) are just the definitions of V, V M and S**" , and therefore are 
in our approach always satisfied. On the other hand Eq. (31) gives in the first order 

(p 2 - M 2 )^!) = 2 [o-(O)0-(l) + 7T(o)7r(i)] .F(o) + d v -K{Q)A\ 0) (78) 

and 

(p 2 - M 2 )V {1) = 2 [o-( O )0-(i) + T(o)T ( i)] P( ) + d v * {0) A\ 0) , (79) 
whereas using Eq. (34) one finds 



(p 2 - M 2 )^ 1} = 2 [a {0) a {1) + 7r (0) 7r (1) ] A[ 0) - M 2 Fd^. (80) 

One can notice that Eq. (80) as it stands is a chirally invariant expression. On the other 
hand Eqs. (78) and (79) form a system of chirally invariant equations in the same way as 
equations (26) and (27), or equations (29) and (30). One can also check, that using formula 
(32) as the definition of V and Eq. (78), we can derive Eq. (79). Therefore, only one of 
equations (78) and (79) is really independent. We observe that all the expressions (78) - 
(80) are just the generalized mass-shell constraints for J~{i), *P(i) an( i *^(i)> an d they do not 
influence the relations in the zeroth order. 

The last equation that we are still required to check is Eq. (38). Substituting the quantity 
V, calculated from (32) into (38), gives 

M 2 ( X) 

—^J-d^(x)d^F(x,p) = rrjg\x,p). (81) 

We see that Eq. (81) defines the parallel part of A^ in the first order. Multiplying (80) by 
p^ and using Eqs. (81) and (54) we obtain zero. Hence Eqs. (81) and (80) are consistent. 
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After deriving all possible equations up to the first order in % and checking their con- 
sistency, we are still faced with two problems. Firstly, the system of our equations in the 
leading order is not closed, since we do not have an equation for Secondly, it is still 

not clear whether our kinetic equation for spin satisfies the requirement of the axial current 
conservation. 

One a priori possible situation is that the condition of the axial current conservation 
eliminates the freedom connected with the choice of d^(X). However, by studying some 
simple situations, one can convince oneself that this is not the case. For example, assuming 
that M = const and that initially the spin distribution is homogenous A^(t = 0,x,p) = 
B^(p) we nn d y^o-^o) = £° l ^ k Pi Bkijp)- In this situation, the axial current conservation 
at t = requires that / d A p £ 0l ^ k (j>j / p°) Bk{jp) = 0. We can see that this equation is 
not sufficient to determine the gradient of $ and consequently the derivative c?o«4° ) remains 
not well defined. Consequently, in order to obtain the closed system of equations in the 
classical limit we have to make some assumption on the form of d^(X). The simplest 
chirally invariant choice is 

d^{X) = 0. (82) 

In the case (82), our kinetic equation (68) allows us to determine the time evolution of the 
function A^(X, p) from its knowledge at some initial time. One can notice, however, that 
this equation allows for solutions which do not satisfy the requirement of the axial current 
conservation. One of these solutions, for the case M = const, has the form A^(X, p) = 
s"(p)s v (p)X v 6^\p - p), where is a vector satisfying the conditions s^(p)p^ = and 

s>1 {p) s h{p) = ~~ 1> an( i P 1S some given value of the four-momentum, with p^p^ = M 2 . On the 
other hand, there also exist solutions of Eq. (68) which satisfy the requirement of the axial 
current conservation. These are, e.g., the static homogeneous solutions (for the case M = 
const) or the solutions which are odd functions of momentum. 

In consequence, we observe that Eq. (68) neither guarantees nor contradicts Eq. (10). 
Therefore the axial current conservation should be used as an external condition which selects 
the physical solutions of (68). 

9. Explicit breaking of chiral symmetry 

In the case m ^ 0, the chiral invariance of the Lagrangian (1) is explicitly broken by 
the expression —m i^i^. In this situation, the additional term in the Dirac equation leads 
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to a simple modification of our kinetic equations (24) - (28); we have to replace the mean 
field a by the sum a + m . Of course, in practice, this change affects only the real parts of 
our equations, i.e., the formulae (31) - (35), since the imaginary parts (37) - (41) depend 
only on the gradients of a. On the other hand, we observe that Eqs. (58) and (59) remain 
unchanged because they are just definitions of the mean fields. 

The important fact in the non-symmetric case is that the relations between the functions 
•7"(o) 5 T'(o) an( i V(p) l°°k different. In particular, instead of Eqs. (52) and (53), we now find 

no) = -m-^—> Vft = p"^2-. (83) 
o-(o) + m ( > o-(o) + m 

Using the first of the relations (83) in Eqs. (58) and (59), we obtain 

m {X) + 8gJ £?-T {0) (X,p) = 0, (84) 



(27T) 

and 



7T ( o)(X) 



o- (0) (X) +m + SG J ^-T {0) (X,p) 



0. (85) 



(27T) 

A straightforward consequence of these two equations is the condition 



T(o)(*) = 0, (86) 

which implies also that V(o(X, p) = and $(X) = 0. Therefore, we observe that if the 
system is not chirally invariant, equations (58) and (59) allow for the determination of the 
angle $(X). 

For m 7^ 0, the mass-shell constraints (54) and (57), as well as the kinetic equations 
(62) and (68) preserve their form, only in this case we should use the definitions 

F(X,p) = ^ , M(X) = a {0) + m . (87) 
0"(o) + mo 

The gap equation (65) has now the form 

M = m„ + 4GM J ^ ' [1 - m , p) - /-(X, p)] . (88) 
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Let us now discuss the behavior of the axial current in the case when m 7^ 0. In this 
situation, we no longer have an axial current conservation law, and using Eq. (32), we can 
write 

The integration of (89) over momentum gives 

MfaPO = ~^nd x )- (90) 

In contrast to the m = 0, case we observe that such an integration does not lead to any 
constraint for itself. Therefore, without any restrictions on d^A^ we can treat Eq. (89) 
as a defining equation for 'P(i). One can still check that Eq. (89) leads to the mass-shell 
constraint for V(i) which agrees with (79). The latter formula, for the case m ^ 0, should 
be rewritten in the form (p 2 — M 2 )V{\) = d„MA^ y We note that Eq. (90) describes the 
partial conservation of the axial current (PCAC) in our case. 

10. Summary and Conclusions 

In this paper, the chiral symmetry concepts have been explicitly included in the con- 
struction of the transport equations for quark matter. Our starting point was a chirally 
invariant Lagrangian and we derived the transport equations via a spinor decomposition 
of the Wigner function and a gradient expansion. Our calculation was restricted to the 
mean field approximation. We have taken into account the spin dynamics and discussed the 
possibility of having a non-zero pseudoscalar condensate. In this aspect, our results are a 
generalization of some earlier results obtained in Ref. [11]. 

The classical quark distribution functions satisfy the kinetic equations of the standard 
form with the effective chirally invariant mass M 2 (X) = tt 2 ^(X) + cr^(X), where tt^(X) 
and CT( )(X) are leading terms in the classical expansion of the pseudoscalar and scalar 
condensates. Furthermore the angle $(X), defined through the relation 7T( )(^)/o"(o)(^) = 
tan$(X), must be a constant. However, its value remains undefined, which reflects the chiral 
symmetry of the problem. 

The classical equation for the spin evolution, which has been derived here for the first 
time, is also invariant under chiral transformations. However, its solutions are constrained by 
the additional condition of the axial current conservation, that is not simple to incorporate. 
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The inclusion of the quark mass term into the Lagrangian explicitly breaks the chiral 
invariance. Consequently, we find that as well as tt^(X) must be zero. Nevertheless, 

the general form of the kinetic equations does not change and we can still use them with 
the substitution M(X) = a^(X) + m , where m is the current quark mass. Moreover, in 
this case the requirement of the axial current conservation reduces to a form familiar from 
PCAC. 

Acknowledgements: We thank M. Volkov, F. Witte and P. Zhuang for clarifying discus- 
sions. 
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